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Oh! 

O . The luminosity determination of electron-positron colliders operating in the region 

of low-lying hadronic resonances {Ecm — 1-10 GeV), such as BEPC/BES, DA<^NE, 
KEKB, PEP-II and VEPP-2M, requires the precision calculation of the Bhabha 
^ . process at large scattering angles. In order to achieve a theoretical accuracy at a 

^ I few 0.1% level, the inclusion of radiative corrections is mandatory. The phenomeno- 

logically relevant effect of QED corrections is taken into account in the framework 
of the Parton Shower (PS) method, which is employed both for cross section cal- 
culation and event generation. To test the reliability of the approach, a benchmark 
calculation, including exact 0{a) corrections and higher-order leading logarithmic 
contributions, is developed as well and compared in detail with the PS predic- 
tions. The effect of 0{a) next-to-leading and higher-order leading corrections is 
investigated in the presence of realistic event selections for the Bhabha process at 
the $ factories. A new Monte Carlo generator for data analysis (BABAYAGA) is 
presented, with an estimated accuracy of 0.5%. Possible developments aiming at 
improving its precision and range of applicability are discussed. 
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1 Introduction 



The accurate determination of the machine luminosity is a fundamental ingre- 
dient for the successful accomplishment of the physics programme of electron- 
positron (e~^e~) colliders operating in the region of the low-lying hadronic reso- 
nances, such as BEPC/BES (Beijing), DA$NE (Frascati), VEPP-2M (Novosi- 
birsk), as well as for the BELLE and BABAR experiments around the T at 
KEKB (KEK) and PEP-II (SLAC). In particular, the precise measurement 
of the hadronic cross section at the $ factories requires a luminosity deter- 
mination with a total relative error better than 1% [1-4]. As well known, 
the luminosity L of e+e~ colliders can be precisely derived via the relation 
L — N/ath, where N and ath are the number of events and the theoretical 
cross section of a given reference reaction, respectively In order to make the 
total (experimental and theoretical) luminosity error as small as possible, the 
cross section ath of the reference process should be large, in order to keep the 
statistical uncertainty small, and calculable with high theoretical accuracy. 

At e+e" machines operating in the energy range 1-10 GeV, the best candidate 
fulfilling the above criteria is the process e+e~ — > e+e~ (Bhabha scattering) 
detected at large scattering angles, say in the angular range 20° < < 160°. 
For example, at the $ factory DA<I>NE the KLOE detector can be used to 
detect such events [4], where the Bhabha scattering cross section is significant, 
being of the order of 10"^ nb at a center of mass (cm.) energy around the $ 
resonance {a/s ~ 1 GeV). 

Therefore, on theoretical side, precision calculations of the large-angle Bhabha 
(LABH) cross section are demanded, with a theoretical accuracy at a few 0.1% 
level. This implies to include in the calculation all the phenomenologically 
relevant radiative corrections, in particular the large effects due to photonic 
radiation. Furthermore, such effects should be implemented and accurately 
simulated in event generators, which are strongly demanded by the experi- 
mental analysis. 

At present, the status of the theoretical predictions and generators of inter- 
est for the LABH process at low-energy e~^e~ machines can be summarized 
as follows. An exact 0{a) generator, based on the calculation of ref. [5] and 
modified to match DA<I>NE characteristics, is used in Monte Carlo (MC) stud- 
ies by the KLOE collaboration [4,6]. An independent 0{a) generator is also 
in use by the CMD-2 and SND experiments at VEPP-2M [7]. In both the 
programs used in such MC simulations the effect of higher-order corrections 
is not taken into account. A semi- analytical calculation of the cross section 
for large-angle QED processes, i.e. e~^e~ — » /i+yU^, e"'"e~, 77, below 3 GeV was 
performed in ref. [8]. It includes exact 0{a) plus leading logarithmic (LL) 
higher-order corrections. This formulation is available as a computer code de- 
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scribed in ref. [9]. The recently developed Bhabha generator BHWIDE [10], 
based on the Yennie-Frautschi-Suura approach for the treatment of QED ra- 
diation, appears in the list of simulation tools presently under consideration 
by the BABAR collaboration at the B factory PEP-II [11]. A QED Parton 
Shower (PS) algorithm, which is employed in the present study, is adopted in 
refs. [12,13] for the computation of radiative corrections to LABH scattering. 
These calculations, however, are optimized to high-energy LABH and differ, 
as it will be discussed, in some aspects from the present implementation of 
the PS model. A complete inventory of existing calculations and programs, 
for both small- and large-angle Bhabha scattering, used at very high-energy 
e+e~ coUiders can be found in ref. [14]. 

The paper is organized as follows. In Sect. 2 the theoretical formulation, which 
is based upon a QED realization of the PS method to account for radiative 
corrections due to photon emission, is described. The steps and kinematics 
of the algorithm are reviewed. Sect. 3 is devoted to the description of a new, 
original PS generator for the simulation of the LABH process and based upon 
the formulation previously discussed. A first sample of numerical results from 
the PS Bhabha generator is also given, with particular emphasis on the sim- 
ulation of the Bhabha process at the $ factories in the presence of realistic 
event selections (ES). In the following Sections tests of the reliability of the 
PS approach are shown and commented, both at the level of integrated cross 
sections and differential distributions. To this end, the calculation of the ex- 
act 0{a) cross section and its matching with higher-order LL corrections is 
addressed in Sect. 4. This is meant as a benchmark calculation, developed in 
order to check the physical+technical precision of the PS generator. Detailed 
comparisons between the PS predictions and the results of the benchmark 
computation are given in Sect. 5. Conclusions, open issues and possible devel- 
opments are discussed in Sect. 6. 



2 Theoretical approach and the Pcirton Shower method 

In order to approach the aimed theoretical accuracy, the calculation of the 
QED corrected Bhabha scattering cross section and the relative event gener- 
ation is performed according to the master formula [15]: 



which is based on the factorization theorems of (universal) infrared and collinear 
singularities. Equation (1) can be worked out within a QED PS algorithm 




(1) 



'cm 
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for the calculation of the electron Structure Function (SF) D{x, Q^), both for 
initial-state radiation (ISR) and final-state radiation (FSR). In eq. (1) dao/dQ 
is the Born-like differential cross section relevant for centre cm. energy be- 
tween 1-10 GeV, including the photonic s- and t-channel diagrams and their 
interference, and the contributions due to exchange of vector resonances, such 
as $, J/^ and T. Following the standard procedure described in ref. [16], the 
contribution of hadronic resonances is taken into account in terms of their 
effective couplings to the electron. At cm. energy around 1 GeV, the total $ 
contribution amounts to ^ 0.1(0.3)% for 20°(50°) <^< 160°(130°), where 
is the electron scattering angle. For higher energies, as in the case of BELLE 
and BABAR experiments around the T, the effect of the T resonance to the 
Bhabha cross section is of the same order. It is worth noticing that for the 
low-energy colliders the LABH cross section is largely dominated by t-channel 
photon exchange; hence, its leading dynamics is quite similar to small-angle 
Bhabha (SABH) at LEPl/SLC and LABH at LEP2 and higher energies [14]. 
In the hard-scattering cross section, the relevant correction due to vacuum po- 
larization is taken into account as well, by adopting the parameterization of 
ref. [17]. In the evaluation of the hadronic contribution to the vacuum polar- 
ization, the Euclidean value of the momentum in the photon propagator has 
been used as the appropriate scale for time-like momenta [18]. The effect of 
the running coupling constant at ^/s = M$ is to enhance the cross section by 
^ 2(2.5)% for 20°(50°) <^ < 160°(130°). The factor J{x_, x+,^iab) in eq. (1) 
accounts for the boost from the cm. to the lab frame due to emission of un- 
balanced ISR from the electron and positron legs, while Q{cuts) represents 
cuts implementation. 

The basic ingredient of eq. (1) is the electron SF D{x,Q'^), which represents 
the probability density of finding "inside" a parent electron an electron with 
momentum fraction x and virtuality Q^. It can be explicitly obtained in QED 
by solving the Dokshitzer-Gribov-Lipatov-AltareUi-Parisi (DGLAP) evolution 
equation [19] in the non-singlet channel: 

X 

where P+ (x) is the regularized e — > e -|- 7 splitting function 

2 1 

P+(x) = - S(l -x) j dtPit). (3) 

1 — X J 


Notice that the energy scale entering the SF will be, in general, dependent 
on the specific process under study. No exact analytical solution of eq. (2) is 
known in the literature. After the efforts undertaken in the last decade for 
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the program of precision physics at LEP/SLC, the theoretical situation can 
be summarized as follows (see for instance ref. [20] and references therein): 



(1) approximate (accurate) analytical solutions in the coUinear limit [21]. 
This kind of SFs is implemented in most of the programs for data analysis 

at LEP/SLC [20]; 

(2) exact numerical solution, obtained via a numerical calculation of the 
Mellin transform of the SF. This can be considered as a benchmark solu- 
tion, but it is unusable from the practical point of view for implementation 
in computational tools; 

(3) exact MC solution, obtained by means of the PS approach, which is par- 
ticularly powerful for exclusive event generation [22-24]. The PS method, 
originally developed and widely applied in perturbative QCD [22], has 
been recently introduced in QED [23,24] as a convenient framework to 
compute photonic radiative corrections in e+e~ collisions. In fact, expo- 
nentiation of soft photons and the contribution of multiple emission of 
hard coUinear photons can be automatically accounted for. 

Let us summarize the basics of the PS method. The starting point of the PS 
approach is the Sudakov form factor [25] : 



n(si,S2) = exp 



a r as 



271 



dzP{z) 







(4) 



which represents the probability that an electron evolves from virtuality —S2 
to virtuality — Si with no emission of photons of energy fraction greater than 
e — 1 — x+, where e is an infrared regulator. In terms of the factor n(si, S2), 
the DGLAP equation can be written in iterative form as: 



D{x,s) = Il{s,m^)5{l-x) 
, ds' 



+ / n(s,s')^n(.',m2)|^ / dyP{y)5{x-y) 



T? 

s 



+ / n(.,/)^ /n(/,.")^n(.",m^)x 

J S J s 



(^) / / ^^2-P(a:^i)-P(a;2)5(x - X1X2) H (5) 



where here = s is understood, s = 4E'^ being the total cm. energy. 
Equation (5) suggests the steps to compute D{x,s) by means of a MC al- 
gorithm [22,23]: 
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X 



Fig. 1. X distribution of the electron SF at = 190 GeV. Solid line: numerical 
solution of DGLAP equation, by means of numerical inversion of Mellin transform. 
Histogram: result of the PS algorithm. 

(1) Set initial values for the electron virtuality and momentum fraction: = 
ml and x — 1 (with Wg electron mass). 

(2) Choose a random number ^ in [0, 1] . 

(3) If C < n(s, K^) then stop: no photon is radiated. 

(4) Else if ^ > n(s,fC^), a photon is emitted: calculate the value K'^ of the 
electron virtuality after the branching as the solution of the equation 

(5) Choose the residual momentum fraction z of the electron in [0, , ac~ 
cording to P{z). 

(6) Replace x by zx and K'^ by K'"^. Go back to step 2. 

In this way the emission of a shower of photons by an electron is simulated, 
and the x distribution of the PS event sample reproduces D{x,Q'^). Such a 
distribution, obtained from a 10^ event sample at = s = (190)^ GcV^, is 
compared in Fig. 1 with a numerical solution of eq. (2), normalized to the same 
number of events. The agreement is excellent. A further test of the algorithm 
is the comparison between the exact analytical Mellin moments D{N) of the 
SF and the corresponding ones calculated in the PS scheme. This comparison 
is shown in Fig. 2, where it can be seen that the PS simulation (markers) 
for = 1,2,10,50,200 well agrees, within the statistical errors, with the 
analytical moments (solid line). The results shown in Fig. 1 and Fig. 2 have 
been obtained with the value e = 10~^ for the infrared cut-off entering eq. (4). 
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Fig. 2. Comparison for the Mellin moments of the electron SF at yfs = 190 GeV. 
Sohd hne: exact analytical moments. Markers: results of the PS algorithm for 
A^ = 1,2,10,50,200. 

However, independence of the PS predictions for the QED corrected cross 
section from e variation, on a wide range of e values (from 10"''^^ to 10"^), has 
been successfully checked, as shown in Fig. 3 for the LABH cross section at the 
$ factory. In the present implementation of the PS model e is taken constant, 
in order to avoid loss of accuracy in the determination of the absolute value 
of the electron SF, as first pointed out in ref. [23]. 

An up to 0{a) PS algorithm has been developed as well. It allows to cal- 
culate the corrected cross section of eq. (1) up to 0{a). Such a calcula- 
tion is strongly required for fully consistent comparisons between the PS 
predictions and an exact perturbative calculation. The steps required for 
the 0{a) PS can be obtained by using eq. (5) and expanding the product 
D{x.,Q^)D{x+,Q^)D{y+,Q^)D{y_,Q^) present in eq. (1) up to 0{a). It is 
easy to see that: 

[D{x^, s)D{x+, s)D{y+, s)D{y_, s)]o(a) = 

[n^(s, m2)]o(„)5(l - x+)S{l - x_)S{l - y+)S{l - y_) 

+ 2^ ^^^^^ ~ ""-^^^^ ~ y-"^^^^ ~ ?/-)^(^+) 
+S{l-x+)S{l-y+)S{l-y_)P{x_) 
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Fig. 3. QED corrected Bhabha cross section as a function of the infrared regulator 
e, at the peak of the $ resonance. Cuts used are given in the text and correspond 
to a reahstic event selection at DA$NE. 



+S{l~x+)5{l-x.)S{l-y+)P{y^) 
+5{1 - x+)S{l - x^)S{l - y-)P{y+)}, 



(6) 



where [Il'^{s,m^)]o(a) is the product of four Sudakov form factors expanded 
up to 0{a). 

The steps for 0{a) PS are therefore the following: 

(1) Set initial values for the fermions virtuality and momentum fractions: 

= ml and x+ = a;_ = y^ = y^ = 1. 

(2) Choose a random number ^ in [0, 1]. 

(3) If ^ < [n^(s, K'^)]o{a) then stop: no photon is radiated. 

(4) Else if ^ > [n^(s, K'^)]o{a): a photon is emitted: calculate the value K'"^ of 
the fermion virtuality after the branching as the solution of the equation 

(5) Choose randomly the fermion which has emitted the photon. 

(6) Choose the residual momentum fraction z of the fermion in [0, 1 — e] , 
according to P{z) and replace or y_ with z, according to the 
particle which has radiated. 

(7) Stop the algorithm. 

The PS algorithm, both in the all order and 0{a) implementation, offers 
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the possibility to go naturally beyond the strictly coUinear treatment of the 
electron evolution, by generating the transverse momentum p± of electrons 
and photons at each branching. To this end, a definite role for the variable 
X must be chosen. In the QED PS model of ref. [23] x is understood as the 
longitudinal momentum fraction of the electron after photon emission. In the 
present PS simulation, x is chosen as the energy fraction of the electron af- 
ter photon emission (hereafter denoted as E scheme), in agreement with the 
known perturbative results for the photon spectrum in QED [26] and pre- 
vious interpretation in perturbative QCD (as, for example, in the paper by 
Marchesini- Webber in ref. [22]). It is worth noticing that the two prescriptions 
are coincident in the coUinear limit. In the E scheme, the kinematics of the 
branching process e(p) — > e'(p') -|- 7(g) can be written as: 



p^{E,0,p,) 
p'^{zE,p^,p'^) 

q^((l-z)E,-p^,q,), (7) 

with built-in energy and transverse momentum conservation. After having 
generated the variables k'^, k'^ and z by the PS algorithm, the on-shell con- 
ditions p^ = /c^, p'^ = k''^, — 0, together with the longitudinal momentum 
conservation, allow to obtain complete event reconstruction as follows: 



2E 

p'^^zE-^1-^^''^'" 



= (1 - z)E 



2E 

ze - k'^ 

2E 

,,2 



pi^(l-z)(ze-kn, (8) 

at first order in k'^/E'^ -C 1, p^/E^ <C 1. An alternative procedure, followed 
in the literature [24], consists in generating a photon p±_ according to the 
leading pole behaviour 1/p ■ q. As a cross-check of the results obtained by 
means of eq. (8), the method of ref. [24] has been also employed in the present 
implementation of the PS model, finding agreement between the procedures 
for the calculation of the QED corrected Bhabha cross section. 



3 Bhabha generator and first sample of phenomenological results 



In the spirit of the PS approach described above, a new MC event generator 
(BABAYAGA) for simulation of the LABH process at e+e" flavour factories 
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Fig. 4. Comparison between the QED corrected LABH cross section and the Born 
one, in the cm. energy range from 1 GeV to 10.5 GeV. SoHd circles: Born approxi- 
mation; sohd triangles: QED corrected cross section. Cuts are given in the text. 

has been developed. It is a generator of unweighted events, giving as output the 

QED corrected cross section and the momenta of the final-state particles. Both 
ISR and FSR are simulated. In principle, an "arbitrary" number of photons 
can be generated, including their p±. In the standard version, BABAYAGA 
records the momenta of electron and positron, and of the most energetic and 
next-to-most energetic photons generated by the electromagnetic shower. The 
possibility of an up to 0{a) calculation of eq. (1) is inchidcd as an option, 
in order to compare it with the exact 0{a) perturbative results (see Sect. 5). 
Also for the 0{a) branch, particles momenta are reconstructed. 

As a first example, in Fig. 4, the PS corrected LABH cross section is compared 
with the Born-like cross section, as a function of the cm. energy. An energy 
threshold of 0.8 ■ -Ebeam for both electron and positron is required, the accep- 
tance cuts are 20° < '&± < 160° and an acoUinearity cut of 5° is imposed. In 
such a situation, the correction due to QED radiation is of the order of 10% 
at the $ factories and of the order of 20% at the B factories, pointing out the 
need of a careful treatment of photonic corrections in simulation tools of the 
LABH process at flavour factories. 

A further sample of phenomenological results obtained by means of the BA- 
BAYAGA program is shown in Figs. 5-7. The parameters and selection cri- 
teria adopted in the analysis are very similar to those considered in previous 
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Fig. 5. Relative effect of ISR and ISR+FSR on the integrated Bhabha cross section 
at the $ factories. Open circles: ISR only; solid circles: ISR+FSR. Cuts used are 
given in the text. 



MC simulations [4,6] and correspond to realistic data taking at DA$NE and 
VEPP-2M, for cm. energy at the $ peak, i.e. = 1.019 GeV. The energy 
cut imposed on the final-state electron and positron is E^^^ = 0.4 GeV; two 
different angular acceptances of 20° < ^± < 160° and 50° < ^± < 130° are 
considered, with (maximum) acoUinearity cut allowed to vary in the range 
Cmax = 5°, 10°, 15°, 20°, 25°. The energy cut refers to the energy of the bare 
electron and positron, corresponding to a so-called bare ES (sec for instance 
ref. [14]), which is not far from realistic due to the presence of magnetic fields 
as for DA$NE and CMD-2 experiment at VEPP-2M. The tight acoUinearity 
cut $,max = 5° is generally introduced in MC simulations in order to single 
out quasi-elastic Bhabha events [4]. However, it is worth noticing that this 
acoUinearity cut, in association with the high energy thresholds on bare elec- 
trons and positrons, defines a rather severe set of constraints, which can be 
expected to emphasize the effects of QED radiative corrections, especially from 
the final-state. This conjecture is confirmed by the numerical results of Fig. 5, 
where the (relative) effect, with respect to the Born cross section, of ISR only 
(open circles) is compared with the whole effect of ISR and FSR (solid cir- 
cles). Actually, it can be seen that the photon radiation produces a lowering 
of the integrated cross section of the order of 10% and that half of this effect 
has to be ascribed to FSR, showing the need of including FSR for realistic 
simulations. In this simulation, the entering the electron SF is fixed to be 
— s as typical virtuality in initial- and final-state shower. The relevance of 
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Fig. 6. Effect of FSR on Bhabha differential distributions at the $ factories: electron 
energy, electron scattering angle, electron transverse momentum and acoUinearity. 
Dashed histograms: ISR only. Solid histograms: ISR+FSR. Cuts are given in the 
text. 

FSR in such ES is further illustrated in Fig. 6 for the electron energy, electron 
scattering angle, electron transverse momentum distribution and acoUinearity 
distribution as well. It can be noticed that the electron energy and p± are sig- 
nificantly affected by FSR, while this is not the case for the electron scattering 
angle and acoUinearity, which are actually largely dominated, as expected, by 
the X distribution of ISR [15]. In the comparison shown in Fig. 6 the numbers 
of generated events for ISR only and ISR+FSR are consistently normalized to 
the same luminosity. 

A further illustration of the potential of BABAYAGA in event generation 
is given in Fig. 7, showing the energy, polar angle and transverse momen- 
tum of the most energetic photon, as well as the missing mass distribution 

of the event, defined as ^ (p+ + — q+ — — k)'^, where p±{q±) are the ini- 
tial(final) electron/positron momenta and k is the momentum of the most 
energetic photon. A minimum energy cut of 5 MeV is imposed as visibility 
criterion. It can be seen that the expected physical features of photon radia- 
tion are correctly reproduced by the PS generator, noticeably the soft peaking 
behaviour of E^^p\ and missing mass distribution as well as the initial- and 
the final-state colhnear peaks, which are clearly visible in the cos'd^ distribu- 
tion. 
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Fig. 7. Distributions of the most-energetic photon in the LABH process at the $ 
factories: energy, polar angle, transverse momentum and missing mass. Cuts are 
given in the text. 

4 Benchmark calculation 



In order to test the (physical+technical) precision of the PS approach and 
corresponding Bhabha generator, an exact 0{q) perturbative calculation has 
been addressed, by computing the up to 0[a) corrected cross section as follows 

c^£L = 4+v(^7 < ^o) + (^h\e^ > ko, cuts). (9) 



In the above equation the 0{a) soft+virtual part <7s^y is obtained by inte- 
gration over the electron angle of the soft -|- virtual differential cross section, 
which is explicitly given by [27] 

da'fltr = ci(7o { 1 + 2 (/3e + Pint) In ko/E + C^} . (10) 



In eqs. (9) and (10), ko stands for a (small) photon energy soft-hard (fictitious) 
separator, i is an index for the photonic contributions to the Bhabha Born 
matrix element {i = |7(s)p,|7(t)p,7(s)-7(t)), f3e = ^a/niln^s/ml) — 1) is the 
leading collinear factor for initial- and final-state radiation, Pint = 2a/TTlii{t/u) 
is the leading angular factor for initial-final state interference. Furthermore, in 
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eq. (10) Cp — Cp{'i9) are the soft+virtual i^-factors for the three QED chan- 
nels, including box/interference finite terms [27]. The hard bremsstrahlung 

contribution cr^'* is included, in the photon phase-space region above the en- 
ergy cut ko, via the e^e^ — e+e~7 matrix element calculated in ref. [28]. Its 
contribution is computed numerically in the presence of experimental cuts by 
means of the MC method with standard importance sampling technique, in 
order to handle collinear and infrared photon singularities. Needless to say, the 
independence of eq. (9) from the fictitious infrared cut-off ko was successfully 
checked with high numerical precision. 

The comparison between the exact 0{a) calculation and the 0{a) predictions 
of the PS generator is of interest because it allows to evaluate the size of the 
0{a) next-to-leading order (NLO) corrections, which are missing in the PS. 
Moreover, this comparison can be a useful guideline to improve the agreement 
between perturbative and LL PS results, for example, by properly choosing 
the virtuality in the electron SF in such a way that 0{a) NLO terms are 
effectively reabsorbed into the LL contributions. It can be noticed, in fact, 
that by choosing the scale as = st/u in the collinear logs ln((5^/m^) 
generated by the PS method, then one has that ln(Q^/m^) — > ln(s/m^) -|- 
\n(t/u) [29]. If this procedure is applied both to ISR and FSR by choosing as 
maximum virtuality of the electromagnetic shower = st/u, it is possible to 
keep under control, besides the large logarithms from ISR and FSR, also the 
leading angular contribution from initial-final state interference. 

In order to estimate the precision of the PS approach with all order corrections, 
higher-order LL terms must be added to the exact 0{a) cross section in the 
benchmark calculation. The general algorithm recently proposed in ref. [30], 
and there applied to the high-precision computation of the SABH cross section, 
can be advocated, by writing the benchmark cross section as (in the so-called 
additive form) [30] 

^ _ _ ^(«) I ^(«) (-1-1) 



where is the all-order LL cross section as given by eq. (1), a^l is the 

up to 0(a) truncation of the LL cross section, a^xact is the exact perturbative 
0(a) cross section of eq. (9). Collinear SFs as given in ref. [21] are used in 
the calculation of the LL cross sections. Adopting such a procedure, exact 
0(a) corrections are simply matched with LL higher-orders in the collinear 
approximation. Therefore, the cross section of eq. (11) includes exact 0(a) 
0(q;"L"), with n > 2, leading corrections. 

A more accurate factorized form, as motivated and discussed in detail in 
ref. [30], can be supplied, by computing the cross section as 
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{a) {a) {a) 

Ma) _ exact ~ CTlL _ '^NL f-,r,\ 
'^NL = = 5 U^J 

where Cj^l is the whole NLO content of eq. (9). In such a way, the bulk of 
the most important second-order NLO corrections, i.e. the 0{a^L) terms, are 
added to the content of eq. (11) [30]|3]. 

The procedure here shortly sketched works efficiently for cross section calcu- 
lation (not unweighted event generation) and allows, by comparison with the 
full PS predictions, to estimate the overall precision of the PS approach. Fur- 
thermore, it is possible to get a measure of 0(q;"L") and 0{a^L) corrections. 
The benchmark calculation is available in the form of MC integrator program 
(LABSPV), which is a suitable modification of the SABSPV code [31] from 
SABH to low-energy LABH process. 



5 Tests of the approach and further numerical results 



As already remarked in the previous section, the detailed comparison between 
the exact 0(a) cross section of eq. (9) and the up to 0(a) expansion of the PS 
results is a valuable tool to establish the physical precision of the PS approach, 
since the registered difference is due to the NLO corrections left over in the 
pure LL predictions of the PS scheme. The relative difference between eq. (9) 
and the up to 0(a) PS cross section is shown in Fig. 8 for the angular accep- 
tances 20° < ^± < 160° and 50° < '&± < 130° and for two different choices of 
the scale in the PS, i.e. for = st/u and = 0.75 • st/u. The relative 
deviations shown in the figure are plotted as functions of the acollinearity cut. 
For the "natural" choice of the scale = st/u which allows to completely 
reproduce the LL structure of the exact 0(a) calculation, the difference is 
an unambiguous evaluation of the NLO corrections. Such contributions, as a 
priori expected, are important in view of the required theoretical precision, 
especially at the acollinearity value ^max = 5°, being of the order of several 
0.1% (see Fig. 8). However, it is worth noticing that a simple variation of the 

scale from = st/u to = 0.75 ■ st/u significantly reduces the differ- 
ence, which goes from more than the 0.5-1% level down to about 0.1-0.3% at 
^max = 5°, depending on the angular set up. In general, with the adjusted scale 

= 0.75 ■ st/u the difference between the exact 0(a) calculation and the 
PS predictions is within 0.5%. This naive example illustrates how, for a given 
set of cuts, the level of agreement can be substantially improved by a simple 

^ Actually, in ref. [30] a more refined treatment of the factorized cross section is 
given, and it is shown that eq. (12) is an approximation of the full treatment at the 
0.1% level, which is sufficient for the present study. 
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Fig. 8. Relative differences for the LABH process at the ^> factories between the 
exact 0{a) cross section (LABSPV) and the up to 0{a) PS one (BABAYAGA), as 
functions of the acolhnearity cut and for two choices of the scale in the electron 
SF. On the left the acceptance region 20° < ??± < 160° is considered, while on the 
right the acceptance region is 50° < 'd± < 130°. Other cuts are specified in the text. 



redefinition of the maximum virtuality of the electromagnetic shower. Going 
beyond this simple procedure would require a merging between perturbative 
calculation and PS scheme, which is beyond the scope of the present work. 

In order to assess the reliability of the PS approach in the presence of higher- 
order corrections, the comparison between the benchmark and PS calculation 
has been extended to the computation of the exact 0{a) plus higher-orders 
cross section. The relative difference between eq. (11) and the full PS cross 
section is shown in Fig. 9, still as functions of the acollinearity cut. According 
to what discussed about the 0{a) comparison, the scale in the PS is, for defi- 
niteness, fixed at — 0.75 • st/u. Relative differences contained within 0.5% 
are still observed, confirming the equivalent implementation of higher-order 
LL corrections in the PS and benchmark calculation. This difference between 
BABAYAGA and LABSPV in the presence of higher-order corrections can 
be considered as an estimate of the physical precision of the "modified" PS 
approach for the cross section calculation for realistic ES at the $ factories. 

In addition to the evaluation of the 0{a) NLO corrections, for an assessment of 
the theoretical precision, it is important to evaluate the impact of the higher- 
order LL contributions. The size of LL 0(q;"L'*), with n > 2, corrections can 
be derived, as already remarked, by comparing in the benchmark calculation 
the results of eq. (9) and eq. (11) or, equivalently, in the PS scheme the full all- 
order predictions with the corresponding up to 0{a) truncation. Furthermore, 
the difference between eq. (11) and eq. (12) in the benchmark calculation is 
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. The same as Fig. 8 for the 0{a) plus higher orders Bhabha cross section at 
factories. Cuts are given in the text. 



able to provide an indication of the size of 0{a'^L) corrections, thus yielding an 
estimate of the most important second-order NLO corrections left over in the 
PS method. Such an analysis of higher-order effects leads to the results shown 
in Fig. 10, where the relative effect of the above higher-order corrections is 
shown for the natural scale = st/u. It can be clearly seen that LL 0{a^L^) 
corrections are unavoidable in view of the expected theoretical accuracy. In 
fact, their contribution is at 1-2% level at ^rnax = 5° and of the order of some 
0.1% for larger acoUinearity cuts. The impact of 0{a^L) is, instead, negligible, 
being well below the 0.1% level, in agreement with the estimate given in ref. [8]. 

Having established that the inclusion of higher-order contributions does not 
alter, as expected, the results of the 0(a) comparison, further 0{a.) tests have 
been performed at the level of exclusive differential distributions. The results 
are illustrated in Figs. 11-13. In these plots, the histograms represent the 
distributions of the events generated by means of the up to 0{a) PS generator 
BABAYAGA with scale = 0.75 ■ st/u, while the markers correspond to 
the predictions of the benchmark calculation LABSPV. A few comments are 
in order here. Since the benchmark calculation is not suited for unweighted 
event generation, the number of events corresponding to the markers have 
been obtained by calculating, as a first step, the integrated cross section over 
each bin and next by multiplying it for a reference luminosity calculated as 
L — Nps/aps, where Nps and aps are the number of events and the cross 
section obtained with the PS generator, respectively. The set of cuts used for 
the analysis of the distributions corresponds to the angular range 20°- 160°. 
Generally, as it can be seen, the level of agreement is quite satisfactory. In 
particular, the electron energy distribution simulated by BABAYAGA, well 
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Fig. 10. Effect of higher-order 0(a"L") (n > 2) and 0{a^L) corrections on the 
integrated Bhabha cross section at the ^ factories. On the left the acceptance re- 
gion 20° < ■d± < 160° is considered, while on the right the acceptance region is 
50° < '&± < 130°. Other cuts are specified in the text. 



agrees with the exact LABSPV calculation, as shown in Fig. 11. Because the 

electron energy distribution is driven in the present implementation of the PS 
method by the x distribution of the SF, the agreement observed between exact 
and PS predictions reinforces, a posteriori, the interpretation of the x variable 
as residual fraction energy as the most natural one in QED PS models. 

Also the PS description of the electron angular variables is in nice agreement 
with the exact calculation, as shown by the electron scattering angle and 
acoUinearity distribution of Figs. 12-13. Some disagreement is seen in the first 
acoUinearity bins, where still acceptable differences at a few per cent level 
are registered. The satisfactory PS predictions for the electron angles can be 
understood since the distributions of these variables are dominated, beyond 
the tree-level approximation, by the effect of the longitudinal boost due to the 
emission of unbalanced electron and positron ISR, which, in turn, is governed 
by the x distribution of the electron SF. 



6 Conclusions and perspectives 



The LABH process is used at present $ and B factories as the reference reac- 
tion to determine the machine luminosity. In order to reach a total accuracy at 
a few 0.1% level, precision calculations of the LABH scattering cross section 
and distributions become more and more urgent. In particular, the relevant 
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Fig. 11. Electron energy distribution for the LABH process at the ^ factories. Mark- 
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Fig. 12. Electron angle distribution for the LABH process at the $ factories. Mark- 
ers: exact 0{a) via LABSPV. Histograms: PS prediction via BABAYAGA with 
scale = 0.75 • st/u 
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Fig. 13. Acollinearity angle distribution for the LABH process at the $ factories. 
Markers: exact 0(a) via LABSPV. Histograms: PS prediction via BABAYAGA 
with scale = 0.75 ■ st/u 

effects due to photon emission have to be kept under control and accurately 
simulated in the computational tools required by the experimental analysis. 

Along this direction, an original calculation of the LABH process has been 
addressed. It is based on a realization of the PS method in QED to account 
for photonic corrections due to ISR, FSR and initial-final-state interference. 
The approach adopted allows the calculation of integrated cross sections as 
well as event generation, including reconstruction of photon p±. A new MC 
event generator (BABAYAGA) has been developed and is available for a full 
experimental simulation of the LABH process at flavour factories. The pro- 
gram has been used to provide several numerical results of phenomenological 
interest, thus showing the potential of BABAYAGA in physics analysis. For 
example, it has been shown, both at the level of cross section and interesting 
distributions, that the effects due to radiation from final-state electrons have 
to be carefully taken into account in the presence of realistic ES at the $ 
factories DA$NE and VEPP-2M, since the impact of FSR can be as large as 
the one due to ISR. 

With the aim of checking the overall (physical-|-technical) precision of the PS 
approach and relative generator, a benchmark calculation has been carried out 
as well. It relies upon the exact 0{a) calculation of the LABH cross section 
supplemented with higher-order LL corrections in the collinear approxima- 
tion. The benchmark computation is also available in the form of computer 
code (LABSPV), which is a MC integrator allowing for precise cross section 
calculations, with an estimated precision at 0.1% level. 
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By comparing the predictions of the PS generator BABAYAGA and of the 
benchmark calculation LABSPV for the up to 0{a) cross section with typical 
cuts, it turns out that the contribution of the 0{a) NLO corrections is, as 
expected, important in the light of the required theoretical precision, being 
at the 0.5-1% level for typical ES at the $ factories. However, it has been 
also shown that the scale entering in the PS LL calculation can be simply 
adjusted in order to agree within 0.5% with the benchmark calculation. An 
energy scale hke Q'^ ^ st/u reveals to be the best choice that effectively 
reabsorbs 0{a) NLO terms into LL PS contributions. This conclusion, which 
holds for the integrated cross section, has been proved to be generally valid 
also for the most interesting differential distributions. 

With the scale choice Q"^ = st/u, the effect of higher-order 0(q;"L") LL cor- 
rections has been evaluated and found to be at the 1-2% level, while the role 
of 0{a^L) corrections is marginal, below 0.1% accuracy. Therefore, one of the 
main conclusions of the present analysis is that, given the size of the radiative 
corrections discussed above, theoretical predictions for the LABH process at 
flavour factories aiming at a few 0.1% precision must include the contribu- 
tion of both 0{a) NLO terms and 0{q:'"'L'"') leading logarithmic contributions. 
From the whole of the present analysis, it turns out that the present physical 
precision of BABAYAGA generator is 0.5% for typical ES at the $ factories. 
If particularly stringent requirements of theoretical accuracy would be in the 
future necessary, the PS approach as presently implemented in BABAYAGA 
should be updated by means of an appropriate merging of the exact 0{a) 
matrix element with the exclusive photon exponentiation realized by the PS 
algorithm. Some proposals addressing such an issue are already available in the 
hterature [32] , in order to obtain sensible QCD phenomenological predictions 
for experiments at the TEVATRON and LHC. 

A second possible development of the present work concerns application of 
the PS scheme to other phenomenological studies of strong interest at e+e" 
flavour factories. In fact, since the PS is a very general method to com- 
pute photonic radiative corrections, the same formulation here applied to the 
LABH process could be employed to evaluate radiative corrections to other 
large-angle interesting QED processes, such as e'^e~ — >• /x^/x" (717), 77 (717), 
and also to hadronic final states, in particular for e'^e~ — > hadrons and 
e~^e^ hadrons + 7, the latter of great interest for an energy scan of the 
hadronic cross section below the nominal cm. energy [33]. 

A further foreseen development is a phenomenological analysis of QED pro- 
cesses at the B factories, along the hues followed in the present study. 

Such possible developments are by now under consideration. 
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